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I.  INTRODUCTION 


In  this  paper  we  consider  the  behavior  of  single  server  queues  In 
which  the  service  time  distribution  S(t)  Is  exponential.  The  1th  Inter¬ 
arrival  time  Is  assumed  to  have  a  distribution  A^(t)  (or  density  a^(t)). 
In  the  most  general  case  all  A^'s  may  be  different.  Here,  however,  we 
consider  the  case  where  the  sequence  of  A^'s  consists  of  a  repeating 
cycle  of  length  n.  In  other  words  the  kth,  n+kth,  2n+kth,  ...  arrivals 
have  the  same  Interarrival  distribution  A^  for  k»l  ...  n.  These  distri¬ 
butions  are  assumed  to  have  rational  Laplace  transforms.  We  refer  to 
this  queueing  system  as  GnjM|l  system.  Such  queueing  systems  arise  In 
computer  networks  using  deterministic  routing. 

The  behaviour  of  a  Gn|M|l  system  may  be  characterized  in  terms  of 
its  virtual  waiting  time  where  the  virtual  waiting  time  at  any  Instant 
Is  the  unfinished  work  at  the  server.  This  Is  a  well  known  approach 
used  in  the  study  of  GI|G|1  systems  [1,  2,  3,  4]. 

virtual  waiting  time  for  a  G|M|1  queue  Is  known  to  have  an 
exponential  steady  state  distribution  [3,  5].  Here  we  prove  that  the 
steady  state  distribution  for  the  virtual  waiting  time  of  a  Gn|M|l 
queue  Is  a  sum  of  n  exponential  terms.  The  approach  taken  In  this  paper 
uses  the  decay  characteristics  of  the  expected  virtual  waiting  time  In 
the  future.  In  section  II  after  Introducing  the  necessary  notation, 
the  equivalence  of  the  expected  virtual  waiting  time  decay  function  and 
the  virtual  waiting  time  distribution  function  are  established.  The 
solution  for  the  steady  state  expected  virtual  waiting  time  decay 
function  for  a  6n|M|l  queuing  system  Is  obtained  In  section  III.  An 
example  for  application  of  the  results  of  a  Gn|M|l  queue  Is  given  In 


Section  IV. 


II.  NOTATION 


The  approach  used  here  Is  based  on  the  results  presented  In  [4]. 

For  the  sake  of  completeness.  In  the  following  we  present  the  notations 
and  results  relevant  to  the  discussion  here. 

Consider  a  server  providing  service  to  arriving  customers 

according  to  a  first- In- first-out  discipline.  Let  the  service  need  of 
a  customer  be  determined  probabilistically  and  defined  by  a  service  time 
distribution  S(t)  (or  a  corresponding  service  time  density  function  s(t)). 

We  are  Interested  In  characterizing  this  queuing  system  In  terms  of  the 
virtual  waiting  time  x(t)  at  time  Instant  t.  Let  the  density  and  the 
distribution  functions  of  x(t)  be  w(x,t)  and  W(x,t),  respectively.  Clearly, 
these  functions  are  defined  for  t  2  0  and  t  *  0. 

In  the  virtual  waiting  time  density  w(x,t),  T  is  the  random  variable 
and  t  Is  a  parameter  of  this  density  function.  For  a  given  value  of  t, 
t  specifies  the  waiting  time  from  t  onwards  and  hence  t  specifies  the 
origin  for  t.  When  we  change  the  value  of  x  to  tj  in  w(x,t)  the  origin 
for  x  Is  changed  with  it  to  reflect  the  waiting  time  corresponding  to 
the  Instant  t. . 

When  an  arrival  occurs  at  time  T  the  waiting  time  density  function 
Immediately  after  the  arrival,  w(x,T+)  may  be  expressed  as 

w(x,T+)  ■  j £  w(o,T“)  s(x-c)<to  -*  (2.1) 

As  a  result  of  the  arrival,  the  waiting  time  undergoes  a  step  change 
In  that  x(T+)  »  x(T")  +  service  time.  Therefore,  the  density  of  the  waiting 
time  after  the  arrival  at  time  T  is  obtained  by  convolving  w(x,T")  with 
the  density  function  of  the  service  time  $(x).  Similarly,  the  distribution 
function  of  the  waiting  time  after  an  arrival  at  time  T  Is  given  as 
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W(t*T+)  •  f  W(o,T  )  s(T-e)do  *•  (2.2) 

J0 

or 

W(t.T*)  •  H(t,r)0$(t)  *♦  (2-»> 

Where  ®  Is  tho  convolution  operator. 

Haring  soon  how  w(r»t)  changes  whan  an  arrival  occurs*  lot  us  soo  how 
It  changes  with  tine  when  no  arrival  occurs.  Given  w(t»T)  and  that  no 
arrivals  occur  between  T  and  T+T^*  w(t*T+Tj)  Is  given  as 

w(t,T+T| )  •  5(t)  J wlatT)do  +  wfr+T,  ,T) 

°  t  »  0  h  (2.4) 

and 

WU-.T+Tj )  •  .T)  t  >  0  -»  (2.5) 

where  Sir)  Is  the  unit  inpulse  function  [3]. 


the  figure.  The  magnitude  of  the  delta  function  at  t»0  for  w(t,  T+Tj) 
corresponds  to  the  probability  of  the  event  (0  st(T) <Tj}.  This  pro¬ 
bability  Is  simply  the  shaded  area  of  w(t,T).  The  corresponding  changes 
to  the  distribution  function  (Equation  2.5)  are  shown  in  Figure  2. 

If  we  are  given  the  arrival  Instants  and  the  Initial  waiting  time 
distribution,  we  may  compute  w(r,t)  for  any  time  t  using  equations  (2.1)- 
(2.5).  Note  that  no  assumptions  are  made  about  the  characteristics  of 
the  arrival  process.  The  complete  transient  solution  for  the  waiting 
time  distribution  may  thus  be  obtained. 

If  we  are  only  Interested  In  the  epochs  just  before  or  just  after 
the  arrivals  and  not  necessarily  at  particular  time  Instants  we  may  use 
the  following  result: 

i.L 

Let  Wn+(T«0)  be  the  waiting  time  density  function  after  the  n 
arrival  at  time  0.  The  waiting  time  density  function  just  before 
the  next  arrival  n+1  Is  given  as 

w  -(t)  *  /*wn(t,t)  a_(t)dt  -(2.6) 

(n+1)  Jo  "  n 

where  an(t)  Is  the  Interarrival  time  density  for  the  next  arrival 
and  wn(t,t)  Is  computed  from  wn(x,0)  using  equation  (2.4). 

These  results  can  also  be_used  to  .ohtain.the  steady  state^solution. 

If  one  exists.  To  arrive  **t  the  steady  state  behavior  we  observe 
that  the  waiting  time  distribution  remains  the  same  Immediately  before 
(or  Immediately  after)  each  arrival.  Let  an  arrival  (say  the  nth)  occur  at 
time  0  and  the  system  Is  In  steady-state.  Then,  given  the  next  arrival 
occurs  at  t,  w(~,t)  Is  obtained  using  equation  (2.4).  But  as  random 
arrivals  may  occur,  w(t,t)  Is  only  a  conditional  density  function  for 
the  waiting  time  seen  by  the  next  arrival.  This  marginal  density  is  given 
by  equation  (2.6).  Under  the  steady  state  condition,  we  may  equate  either 


w^Ct)  Wn+1"^T^*  or  wn+^  40  wn+l+^*  Therefore»  W(T)  *s  the  steady 
state  waiting  time  distribution  function  If  and  only  If  W(t)  satisfies 

/*t  rm 

W(t)  ■  /  /  W(c+t)  a(t)  s(T-o)dtdo  -*  (2.7) 

JQ  J0 

Note  that  by  rearranging  the  terms  In  equation  (2.7)  Llndley's  Integral 
equation  [2]  results. 

While  x(t)  Is  a  stochastic  process  which  is  completely  characterized 
by  the  deslty  function  w(x,t),  In  many  practical  situations  we  are  interested 
In  calculating  the  expected  value  of  x(.t),  E^Ct).  Knowing  w(x,t),  E^t) 
may  be  calculated  using  the  following  equation 


Ew(t)  t  w(f,t)dt 


(2.8) 


Note  that  equation  (2.8)  Involves  w(x,t).  In  order  to  use  it  in  calculating 
E^(t)  for  an  arbitrary  t,  w(x,t)  must  be  known  for  all  x. 

We  define  E^tt)  and  E^U),  to  be  the  expected  virtual  waiting  time 
at  Instant  t  evaluated  immediately  after  and  before  an  arrival  occurlng  at 
t*0,  respectively. 

As  shown  in  [4],  E^t)  before  and  after  an  arrival  Is  related  as  follows. 
E>(+(t)  -  Es(t)  +  Ew.C0)Cl-S(t))  ♦  s(t)  ©  Ew.(t) 

-  (2.9) 
and 

Ew_(t)  -  /"  a(T)  Ew+(t+T)dT  t  a  0  ->  (2.10) 


where 


Es(t) 


L 


T  s(t+T)dT.  -»  (2.11) 

Equation  (2.9)  allows  us  to  compute  E^ft)  given  E^_(t)  and  the  occurence 
of  an  arrival.  Equation  (2.10)  may  then  be  used  to  evaluate  ^.(t)  just 
before  the  next  arrival  given  the  Interarrival  density  function  a(t)  for 
the  next  arrival.  These  equations  can  be  used  to  evaluate  the  steady  state 


values  for  these  functions  when  such  a  steady  state  exists. 


The  relationship  between  Eg(t)  and  w(t,0)  Is  established  In  Theorem  1. 
Theorem  1  -  E^t)  and  w(t,0)  are  equivalent  In  that  given  one  the  other 
is  uniquely  determined. 

Proof  -  (a)  Given  w(t,0),  equations  (2.8)  and  (2.5)  uniquely  determine 


yti. 

(b)  Given  E^t)  we  may  write  [5] 

V'J-  f  [1-H(f ,t)]dT. 

From  equation  (.2.5) 

Ew(t)  *  j*  [l-W(t+r,  0)]dr 


*  f  [1-W(v, 


0)]dv. 


Let  1-W(v,0)  »  U(v,0) 

Thus,  E^t)  *  /  U(v,0)dv. 

Also,  EJt+At)  *  f  U(v,0)dv. 

J  t+At 


Therefore , 


Ew(t+At)  -  Ew(t) 
At 


1 

At 


t+At 

U(v,0)dv 


In  the  limit  as  At.  -*  0 


dE  (t) 

— “ —  -  -  U ( t , 0 )  -  -(1-W(t,0)). 
dt 

dEjt) 

Thus  W(t,0)  «  — * —  +  1. 
dt 

This  completes  the  proof  of  'Theorem  I.  ■ 

The  approach  taken  In  this  paper  uses  Ew(t)  to  characterize  the  queuing 
system.  Clearly,  from  E^(t)  the  distribution  functions  for  the  virtual 
waiting  time  can  be  evaluated. 


III.  SOLUTION ’TO  Gn|M 1 1  SYSTEM 


Consider  a  flrst-ln-first-out  exponential  server  with 
s(t)  »  ue’ut. 

Let  the  1“*  Interarrival  distribution  and  density  be  A^(t)  and  a^(t), 

respectively.  Also,  let  E^_  (t)  and (t)  be  the  expected  virtual  waiting 

^"1  '  **1 

time  decay  functions  just  before  and  after  the  1  arrival. 

Assume  that  the  arrivals  to  this  server  consist  of  a  repeating  cycle 


of  length  n  In  that 

*1+kn(t>  "  ai(t> 

In  the  steady  state.  If  one  exists, 

E*+  (t)  *  (t) 

^  1+kn  ^1 


Ew_  (0  »  (t). 

"  1+kn  w  1 


for  1  *  1,2,. ..,n 
and  k  *  1,2,3 _ 


for  1  *  1,2,. ..,n 
and  k  *  1,2,3.... 


-  (3.1) 


To  solve  this  system  we  use  equation  (2.9)  to  evaluate  (t)  given 

i  (t)  and  use  equation  (2.10)  to  evaluate  (t)  given  (t). 

w-i  "“1+1  w+i 

These  combined  with  equation  (3.1)  are  then  used  to  obtain  the  steady 

state  values  of  Ew( t)  functions.  Note  that  for  the  exponential  service 


Es(t)  •  i  e"ut 


S(t)  -  1  -  e'ut 


From  equation  (2.9) 


Also, 


E*.  (t)  ■  ie'ut  ♦  E^.  (0)(e_ut)  ♦  Ew_  (t)  ©  u< 

V*/‘)  '  'Ew-,'0)  *7>*'“*  *  ®  “‘'Ut 

t 

Ew-i(0)  +  v  “  W°> 


t  *  0 


Thus, 


Ew+iU)  -  CO)e"yt  +  Ew_^(t)  ©  ue'ut  -  (3.2) 

t  i  0 


From  equation  (2.10) 

E w-i+1(t)m/0  ai+l(T)  Ew+1(t+T)dT  ta0  -(3.3) 


Let  us  define  the  Laplace  transforms  of  Ew.  (t)  &  E^  (t)  as 

w+i  W-.j 

^  OQ 

Lw+i(s)  =  e”stEw+^  (t)dt  s  a  0  -»  (3.4) 

and  L  (s)  *  f  "  e'SV  (t)dt  s  *  0  -  (3.5) 

w  i  J  0  w'i 

Note  that  as  E.  (t)  and  E  (t)  are  defined  for  t  iO,  L.  (s)  and  L.  (s) 
w  1  w”i  w  i  w"i 

are  analytic  in  the  right  half  plane  of  s.  Also,  as  E  (t)  functions  are 

w 

positive  and  integrable  the  transforms  L  (s)  exist  [3]. 

Talcing  the  transform  of  equation  (3.2)  we  may  write 

Lw+^s^  =  I!+s  +  Lw-^s^  '  ■£+?  *  ^3*6^ 

Equation  (3.3)  expresses  a  form  of  convolution  which  is  defined  only 
for  t  i  0.  If  we  were  to  add  a  function  of  the  form 


f‘  •i+i<T>e1*1<fTWT 

we  may  wri te  the  correspondi ng  tranf orms  as  [ 3 ] 


t  <  0 


G"i+i(s)  +  L  (s)  *  L  (-s)  •  L  (s) 
i+i  w-i+1  ai+1  w+i 


(3.7) 


where  G”^+^(s)  is  a  function  of  s  which  is  analytic  in  the  left  half  plane 
of  s. 


Before  proceeding  further  we  need  the  following  Lemmas. 

Let  a(t)  be  the  density  function  of  the  Interarrival  time  with  a 
rational  laplace  transform  L_ (s ) . 

a 

Lemma  1  Lj-s)  can  be  written  as 
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V, 

Ji 

1 


(3.8) 


La(-s)  *  *+(s)+V(s) 

where 

*+(s)  Is  analytic  for  Re(s)  >  0  and  $"(s)  is  analytic  for  Re(s)  <  0  . 
And  that  $+(s)  has  the  form  $+(s)  *  ~  where  c  is  a  constant. 

Proof.  As  L.(s)  is  rational  we  may  write  it  as 

_  a 


n-1  ,  _n-t  ,  , 

L  (s)  =  VIs  +  VftS  +  —  *ao 

a 

0ns"  +  6n-lsn  1  +  +  s0 

Note  that  a(t)  is  defined  for  t  >  0  only.  Thus,  La(sJ  is  analytic  in  the  right 
half  plane  and  has  poles  onTy  in  the  left  half  plane. 

Therefore,  L  (-s)  is  analytic  in  the  left  half  plane  and  has  pole*  only  in  the 
right  half  plane.  Clearly,  the  product 

has  one  pole  in  the  left  half  plane  and  the  rest  in  the  right  half  plane. 

Therefore,  this  function  can  be  written  as  <5,+(s)  +  $  (s)  where 

$+(s)  *  -TT  where  c  is  a  constant, 

u+s 

This  completes  the  proof.  ■ 


Lemma  2 


where 


(\n  n  +  _ 

-H-  )  n  L.  (-s)  can  be  written  as  $  (s)  +  *  (s) 
u+s  /  1»1  a1 


n  c 
«+(s)  -  I 


1 


T 


(3.9) 


1-1  (u+s) 

Proof.  This  directly  follows  from  the  Lemma  1.  ■ 

Me  observe  that  multiplying  the  expression  by  an  rational  function 
which  is  analytic  in  the  left  half  plane  will  still  yield  a  decomposition 


In  which  ♦+(s)  has  the  form  shown  In  equation  (3.9) 


Lemma  3 


where 


L*.n(s)  ■  x(s)Lw_^(s)  +  Y(s)  -  z(s) 


n  i  I  n  La/'s^ 

Y(s)  -  L  (0)  {  n  -4- - 

1-1  ^1-1  (  j-1  y+s 

n-1  (  n  u 

z(s)  •  I  j  n  -H_l  (-s)  G'(s) 

1»1  ( j-1+1  w+s  aj  J  1 


(3.10) 


Proof.  Follows  from  a  repetitive  application  of  equations  (3.6)  and  (3.7), 

For  the  next  two  Lemmas  let  us  consider  a  distribution  function  Fix'), 
x  as  0  with  expected  value  S  arid  Laplace  transform  JPf(s) . 

Lemma  4  If  (a)  Re(s)  as  0,  |w|<  1  or  (b)  Re(s)  >  0,  |w|s  1  or 

(c)  u6  >  n  and  Re(s)  *  0,  |w|  je  1  then  the  equation 

zn  *  vfcCf(s  +  v(l-z)) 

has  exactly  m  roots  z  -  «r(s,w)(r»li2,...m)  In  the  unit  circle  |z|  <  1. 

We  have 

-  J,m\S 


Ms.w)  »  I  - rr-!- - 

j-1 


(3.11a) 


where  ^  ■  e211^"  (r-l,2,...m)  are  the  mth  roots  of  unity. 

Proof.  See  [6,  page  126].  ■ 

?  Let  6r  •  6f(0,l)  where  6f(s,w)  Is  given  by  equation  (3.11a). 

If  u8  >  m  then  *<52*  •  *6m  are  th®  m  roots  1n  z  of  the  equation 


zm  -^f(wO*))  -  (3.11b) 


In  the  unit  circle  |z|  <1.  If  u3  s  m  then  £^>£2,**,5mol  ar*  ***•  m”3 
roots  in  z  of  (3.11b)  in  the  unltclrcle  |z|  <  1,  whereas  6m  »  1. 

Proof.  See  [6,  page  126].  ■ 

We  further  observe  that  equation  (3.11b)  always  has  one  root  at  z  ■  1.  This 
Is  because  forQ^(s)  to  be  a  Laplace  transform  of  a  probability  distribution 


f 


Theorem  2  -  In  a  Gn|M|l  queue  where  the  arrival  distributions  have  rational 
Laplace  transform, 

L.  (s)  -  I  -  (3.12) 

w  1»1  sn1 

where  Lw_(s)  Is  the  Laplace  transform  of  Ew_(t)  at  the  beginning  of  a 
cycle,  a^'s  are  constants  and  c^'s  are  the  n  roots  of  the  equation 


d  -  n  -t-La  (-.)]  -  o 

i»i  M  s  ai 

In  the  region  |s|  <  p. 

Proof  -  For  a  Gn|M|l  queue  with  a  cycle  of  length  n,  In  the  steady  state. 

If  one  exists,  Lw_  (s)  »  L^_  (s)  •  Lw_(s).  In  order  to  find  L^ts) 
n  0 

we  use  Lenina  3  and  the  technique  of  spectral  factorization  [7].  We  may 
write  equation  (3.10)  as 

Lw_(s)  •  x(s)  Lw_(s)  +  Y(s)  “  *(s) 

or  (1  -  x(s))>  Lw.(s)  ■  Y(s)  -  z(s)  -»  (3.13) 

In  order  to  use  the  spectral  factorization  technique  we  need  to  express 
1  -  x(s)  ■ 

*  (s)  -  (3.14) 

where  (s)  Is  analytic  In  the  right  half  plane  and  <p"(s)  Is  analytic  In 
the  left  half  plane.  From  equation  (3.10) 


(s>  ■(&)" 


•*  (3.14a) 


As  L,  Cs )  Is  the  Laplace  transform  of  a  density  function  a,(t)  and  II  L  (s) 
a1  1  1«1  a1 
Is  the  Laplace  transform  of  the  sum  of  n  random  variables  coming  from 

n 

densities  a.(t),  a9(t)  ...  a_(t).  Thus,  n  L,  (s)  Is  a  Laplace  transform 
i  £  n  w  a1 

of  a  valid  density  function.  By  substituting 

“r-1 

we  may  express  the  equation 

1  -  x(s)  *  0  -*  (3.14b) 

In  the  form  of  (3.11b).  From  Lenina  5  this  equation  has  n  roots.  Let 

the  negative  of  these  roots  be  Cj,  1  ■  1,2,..., n.  We  note  that  s  ■  0  Is  also 

a  root  of  this  equation.  Using  these  roots  Me  may  Mrlte  ?+(s)  and  ?*(s)  as 

n 

s  n  (s+tj) 

m+(s)  ■  — •  -n —  -*  (3.15) 


and 


<p'(s) 


(s+u)n 

n 

s  n  (s+;, ) 
i»l  1 

(s+p)n(l“X(s) ) 


(3.16) 


Clearly,  y($)  and  qf(s)  satisfy  the  necessary  conditions  for  spectral 
factorlztlon  [7].  From  equations  (3.13)  and  (3.14)  m e  Mrlte 


From  Lemma  2 


Mhere 


L^(s)  «Ts)  ■  <p"($)  [Y(s)  -  z(s)] 
qT(s)  [Y(s)  -  z(s)]  ■  4+(s)  +  $“(s) 
<P+(s) 


n 

E 


T 


(3.17) 


(3.18) 


4  (3.19) 


1-1  (u+s) 

From  equations  (3.17)  and  (3.18) 

L*.<S),  (S)  -  As)  •  ♦■(») 

A  stralghtfonvard  use  of  the  analytic  continuation  arguments  and  Llouvllle's 
Theorem  [7]  Implies 

l*. (•)**(•)  -  *+(s)  -  4*(s)  -  K 
Mhere  K  Is  a  constant. 
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Therefore, 


Lw-<s> 


K  -  0*(s) 
<p+(s) 


Substituting  from  equations  (3.15)  and  (3.19),  we  have 


n 


(sm)r 


s  n  (s+Cj ) 
1*1  1 


[  K  +  z  — Xrl  -»  (3.20) 

L  1*1  (y+s) ' J 


To  evaluate  K  we  note  that  lim  E  (t)  *  0 

t  -»  *  w" 

Hence, 


11m  SL  (s)  -  0  ■  11m  h  *  E  — Xj-1 

i.-0  w'  s-0  i  L  1-1  (u+s)1  J 


n  (s+c*) 
i-l  1 


This  implies  that 


n  c 


i 


K  *  -  Z  -r 
1»1  u1 

Using  equations  (3.21)  in  (3.20)  and  simplifying  we  get 


(3.21) 


l^ls) 


Vl5"'1  *  dn-2s"'2  *  "•  +  d0 


II 

n  (s+;,) 
i-l  1 


where  d. 's  are  constants. 


1 


Or, 


U»> 


,"s,  (  ^t) 


where  a^'s  are  constants. 

This  completes  the  proof.  • 

Theorem  3  -  A  6n|M|l  queue  has  a  steady  state  solution  If  >  n  where 

n  r « 

8  »  E  f  t  a,(t)dt 
i-l  Jo  1 

Proof  -  From  Lemma  5  we  observe  that  when  yg  s  n  one  of  the  n  roots  of  the 
equation  (3.11b)  becomes  6n  ■  1.  This  root  corresponds  to  tn  •  0.  Thus 
from  Theorem  2  lw-(s)  has  a  term  of  the  form  j  . 


As  a  consequence  11m  s  l  (s),  and  hence  11m  E  _(t)  Is  finite.  This  cannot 
be  true  for  a  steady  state  solution  of  the  queue.  Therefore  the  queue 
reaches  a  steady  state  If  and  only  If  u6  >  n  and  hence  the  theorem.  ■ 

The  existence  of  a  steady  state  has  been  established  In  Theorem  3 
under  the  condition  uB  >  n.  As  6  Is  simply  the  expected  length  of  time 
for  n  arrivals,  this  condition  Is  no  different  from  the  standard  steady 
state  condition  of  p  <1. 

A  G | M| 1  queue  corresponds  to  a  Gn i M 1 1  queue  with  n«l ,  since  the  solutions 
given  above  for  the  Gn|M|1  queue  correspond  to  the  well  known  solutions  of 
G|M|1  queue  when  we  take  n  »  1. 

The  form  of  Lw_(s)  given  by  equation  (3.12)  implies  that  In  steady 
state  Ew_(t)  can  be  expressed  as  a  sum  of  n  exponential  terms  where  the 
exponent  coefficients  are  determined  from  the  roots  of  the  equation 
(3.14b).  From  the  form  of  x(s)  as  given  In  equation  (3.14a),  the  roots 
of  equation  (3.14b)  do  not  depend  on  the  order  of  a^'s  In  the  sequence. 

The  Initial  value,  Ew_(0)  does  depend  on  the  sequence,  however. 

To  obtain  the  steady  state  solution  of  a  Gn|M|l  system,  one  approach 
Is  to  apply  the  spectral  factorization  technique  used  In  the  proof  of 
Theorem  2  to  determine  the  coefficients  a^'s.  The  coefficients  ^ ' s 
are  obtained  as  roots  of  equation (3. 14a) using  the  results  of  Lemma  5. 

An  alternate  approach  to  obtaining  the  coefficients  a/s  Is  to  assume 


and  substitute  this  In  equations  (3.2)  and  (3.3)  .  Also, 

L.  (0)  -  £  ot j 

^  1  j-1  J 

t  n 

and  p  (0)  -  i  +  e 

^+1  M  j.l  i 

Thus  from  equation  (3.2)  we  have 
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Vl) '  ft  *  ji  “j]  ,'yt  *  u,'ut  ®[  j,  VSjt] 


,/i  * ;  asi!,-*** ;  ml.-v 

'»  j-i  'ruJ  j-i  y-;j 


(3.23) 


Using  equation  (3.3)  along  with  equation  (3.23)  we  get 


(t)  “{ u  +  *  F^}e"yt'La  (y) 
W  2  fw  j»l  5j-w/  a2 


j, 


(3.24) 


As  the  form  of  E^(t)  remains  the  same  for  each  step  in  the  cycle,  In 
equation  (3.24)  we  must  have 


i. ;  i£i.0 

“  j-i  vu 


Similarly, 


S(t”  jiajfe)L*2"j,,'v 

E-3(,)  *  j,  S(!J)L*3t'j,*'V 


(3.25) 


(3.26) 


)  •  0 


(3.27) 


After  the  n*  step  we  have 


E*_  M  •  E^.  (t)  •  E  Oj  (^7-)  I  n  L  U,)  jt’V 

^1  ^  n+1  j-1  JVu  ;J/  Ck-1  ak  J  ) 


(3.28; 


yielding  n  linear  equations  of  the  type  (3.27). 

Clearly,  the  equality  of  Ew  (t)  and  E  (t)  holds  where 

w'l  w"n+l 


(=%)  fii  vcW 


for  j  ■  1 ,2. ..n 


(3.29) 


The  t^’s  obtained  as  negative  roots  of  equation  (3.14b)  satisfy  equation 
(3.29).  The  values  of  a^'s  are  obtained  by  solving  n  linear  equations  of 
the  form  (3.25),  (3.27),  etc.  A  solution  fOr  a  Gn | H j 1  queue  may,  therefore, 
be  obtained  by  finding  n  roots  of  equation  (3.14b)  and  solving  a  set  of 
n  linear  equations.  This  approach  Is  simpler  than  the  spectral  factorization 
technique.  In  the  next  section  we  present  an  example  of  obtaining  the 
steady  state  solution  for  a  Gn|M|1  queue. 


-16- 


IV.  APPLICATION 


A  situation  In  which  G^/m/|  queues  occur  Is  when  an  arrival  stream 
Is  split  among  k  servers  using  a  deterministic  sequence,  the  arrivals 
seen  by  each  server  form  a  cyclic  pattern.  For  example  consider  the  2 
server  system  shown  In  Figure  3.  The  arrivals  at  box  0  are 


1  1  2 

Let  the  arrivals  to  D  be  Independent  with  a  general  distribution  F(t). 

Let  us  denote  the  1-fold  convolution,  F^(t)  by 

1  times 

F(i)(t)  -  F(t)  <S>  F(t7  ©77.  ©  F(t) 

The  arrivals  seen  by  server  1  follow  a  cyclic  pattern 
F(1)(t),  F(2)(t), 

Similarly ,  the  arrivals  at  servers  2  follow  general  arrivals  from 

F<3>(t>, 

2 

Queue  1  Is  of  the  G  |  M|1  type  and  can  be  analysed  using  the  results  of 
this  paper. 

In  particular,  let  us  consider  the  arrivals  to  0  be  Polsslon  with 
rate  X ,  l.e., 

F(t)  «  1  -  e*xt 


Let  the  service  rates  of  the  two  servers  be  Uj  end  u2>  For  server  1, 

F<|  (t)  Is  the  same  as  F(t)  and  F2(t)  Is  given  by 
F(t)  ®  F(t) 

Thus,  La^ (s)  and  L^(s)  are  ^  and(x+s)  ’  respectively.  Equation  (3.14b) 
reduces  to 

(5 £?)2(&)3*’  -(41> 

Also,  we  have  the  two  linear  equations 


and 


1 

/l*l 

.  ^2**2 

a  n 

yl 

C,-M, 

(;2-i*iT 

*  u 

1 

yi 

(x  \ 

yl 

V“l 

u,-c, 

U<*. 


2  2 


(c2-U|)  2 


-C-)  \x+c2/ 


(4.2) 
2 


(4.3) 


As  a  numerical  example  let  »  0-2  and  X  *  0*1 . 
In  this  case  we  get 


t1  »  0*1491 3 
■2 


U  »  0*23295 


otj  *  1*41565 
a2  •  -0*10789 

Thus  Ew_(t)  -  1*41565  e‘°*14913t  -  0.1079  e'°-23295t 


From  this  we  note  that, 

6^(0)  -  1  -30776 

Ew+(0)  •  6-30776 

The  before  and  after  arrival  values  of  the  expected  virtual  waiting  time, 
after  the  next  arrival,  are  2-4556  and  7-4556,  respectively. 
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A  solution  for  this  example  obtained  by  an  Iterative  solution  of 
the  transient  equations  has  been  reported  In  [4].  The  solution  obtained 
here  Is  exactly  the  same  as  the  one  reported  In  [4]. 

for  this  example,  the  Interarrival  times  seen  by  server  2  have  the 
same  three  stage  hypo-exponential  distribution.  The  expected  waiting 
time  seen  by  an  arrival  routed  to  this  server  by  the  above  sequence  can 
be  obtained  by  using  the  standard  solution  to  G|M(1  queue.  For  u  *  0.1, 
the  expected  waiting  time  Is  11.9403.  Thus,  the  expected  waiting  time 
for  an  arrival  to  these  two  servers  using  the  deterministic  routing  Is 
8.5679.  If  the  routing  were  done  at  random  with  probability  2/3  for 
server  1  and  1/3  for  server  2  the  expected  waiting  time  would  have  been 
10.00.  It  Is  Interesting  to  note  that  In  this  case  a  deterministic 
decision  rule  does  better  than  a  probabllltlc  rule. 


V. 


SUMMARY 


In  this  paper  the  steady  state  behavior  of  an  exponential,  flrst-ln- 
flrst-out  queue  with  cyclic  arrivals  has  been  analyzed.  This  arrival 
structure  Is  a  generalization  of  G|M|1  queues.  The  well  known  result  of 
G|M|1  queues  that  the  waiting  time  Is  exponentially  distributed  has  been 
generalized  to  the  Gn|Mj1  queue.  This  generalization  yields  an  easier 
technique  for  solving  a  given  Gn|M|l  queue.  An  example  to  show  the  use 
of  the  results  has  been  presented,  as  well. 
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